A New Solar System Population of WIMP Dark Matter 



Thibault Damour^ and Lawrence M. Krauss^ 
^Institut des Hautes Etudes Scientifiques, 91440 Bures-sur-Yvette, France 
and DARC, Obaervatoire de Paris-CNRS, F-92195 Meudon, France 
^Departments of Physics and Astronomy Case Western Reserve University Cleveland, OH 44^06-7079 



00 



Perturbations due to the planets combined with the non-Coulomb nature of the gravitational potential in the Sun imply that 
WIMPs that are gravitationally captured by scattering in surface layers of the Sun can evolve into orbits that no longer intersect 
the Sun. For orbits with a semi-major axis < 1/2 of Jupiter's orbit, such WIMPs can persist in the solar system for > 10® years, 
leading to a previously unanticipated population intersecting Earth's orbit. For WIMPs detectable in the next generation of 
detectors, this population can provide a complementary signal, in the keV range, to that of galactic halo dark matter. 



o 

Q 

o\ ■ 
CNi : 

m ' 
> 

\o ■ 

o : 
oo . 

ON ■ 

i!^ ■ 
^: 
6 ■ 



The effort to uncover the nature of the dark mat- 
ter that dominates the gravitational potential well in 
our galaxy, and almost all known galaxies, is perhaps 
one of the most important ongoing experimental pro- 
grams in cosmology. Among the favored candidates for 
non-baryonic dark matter are so called Weakly Interact- 
ing Massive Particles (WIMPs). Such WIMPs naturally 
arise in most supersymmetric extensions of the Standard 
Model in the form of the lightest supersymmetric part- 
ner of normal matter, the neutralino (e.g. A num- 
ber of ongoing experiments, involving underground low 
background detectors, capable of detecting the energy 
deposited due to the elastic scattering of neutralinos off 
atoms, are underway. The tip of the allowed parameter 
space of neutralinos is just beginning to be probed. 

If a positive detection is made in the next generation of 
underground detectors, it will be important to search for 
signatures that definitively distinguish such a signal from 
the possible background noise due to radioactive contam- 
inants in and around the detectors. Possibilities include 
searches for a potential annual modulation and/or an- 
gular anisotropy of the signal (2|j|,|l|] , or indirect searches 
involving the annihilation products of the neutralino pop- 
ulations captured by the Sun |J|,|,|ll or the Earth §,1. 

The fact that WIMPs can be captured by the Sun 
and planets motivates examining whether various compli- 
cated dynamical histories in the solar system might affect 
the local WIMP density near the Earth (e.g. As we 

describe here, a careful consideration of the WIMPs that 
are scattered by a nucleus in the Sun into gravitationally 
bound orbits indicates that a small population can elude 
subsequent scattering in the solar interior (which other- 
wise leads to concentration in the Solar core, followed 
by annihilation) . Because of the non-Coulomb nature of 
the gravitational potential in the Sun, this population 
involves only WIMPs that scatter near the surface of the 
Sun. These WIMPs can, due to the perturbing influence 
of the planets, diffuse out of the Sun and build up over 
time to produce a density in the region of the Earth that 



may be comparable to the background WIMP halo den- 
sity. The range of scattering cross sections and masses 
for such WIMPs is precisely that which is associated with 
the range of parameter space just below current experi- 
mental limits. Thus, if WIMPs are to be discovered in 
the next generation of experiments, then the population 
we describe here should be significant. 

We shall focus on the sub-population of WIMPs that 
scatter on a nucleus located near the surface of the Sun, 
and thereby lose just enough energy to stay in Earth- 
crossing orbits. These will be susceptible to small grav- 
itational perturbations by the planets. We are inter- 
ested then in the differential capture rate, per energy, 
and per angular momentum, of WIMPs in the Sun, and 
in particular only in the fraction of WIMPs that have 
angular momenta in a small range [Jmim Js] where J5 
is the angular momentum for a WIMP exactly grazing 
the Sun. The standard low-energy differential cross sec- 
tion of WIMPs on nuclei of atomic number A is given 
hy dffA = a°AFl{Q) where Q = i^boforo - ^^aftor 

is the energy transferred during the scattering, Fa{Q) 
is a form factor, and dilcm is the scattering solid angle 
element in the center of mass frame. Note that is 
by definition independent of the scattering angle, and 
takes a value that depends on the details of the WIMP 
particle physics parameters. Generalizing the standard 
Q capture calculation (see 0), one can derive the rate 
with which WIMPs scatter on nuclei with atomic num- 
ber A to end up into bound solar orbits with semi-major 
axis between [a,a-\-da\ (corresponding to [a, a-l-da] with 
a = GArmgun/a), and with specific angular momentum 
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Here, the minimum radius r,„in (perihelion) is defined 
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in terms of the minimum angular momentum J, 
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cape velocity at Tmin, nA,x is the density of atomic tar- 
gets A, and WIMPs X, respectively, Vo is the rms circular 
velocity of WIMPs in the Galaxy, and the "capture func- 
tion" Ka {r, a) involves an integral over the WIMP local 
phase space distribution at the Sun, weighted over the 
scattering form factor . While the general form (used 
in deriving our results) is not particularly illuminating, 
in the limit in which the nuclear form factor is neglected 
in the scattering cross section, and in which we also ne- 
glect the relative motion of the Sun with respect to the 
galactic halo, KA{r,a) takes the following simple form: 
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where = /"j^ . 

To obtain the total rate, a sum over all the (sig- 
nificant) values of A present in the Sun must be ul- 
timately performed. We shall be interested in values 
Jmin — Js = RsVcsciRs), whcrc Rs is the Sun's ra- 
dius, and a ~ 1 AU, i.e. a ^ Gmsun/ (1 AU) ~ v'^ where 
VE = 29.8kni/s is the Earth orbital velocity. Note that 
for such values of a, typically a ^ '^csc^ ^^'^ from 

the above formula it is clear that the function _ftr^(r, a) 
is nearly independent of a. 

The Sun scattering events create a population of solar- 
system bound WIMPs, moving (for a ^ 1 AU) on very 
elliptic orbits that traverse the Sun over and over again. 
For the values of WIMP-nuclei cross sections we shall 
be mostly interested in here (corresponding to effective 
WIMP-proton cross sections in the range 4 x lO""*^ — 4 x 
lO^^'i cm^), the mean opacity of the Sun for orbits with 
small perihelions is in the range 10"'' -10^3. This means 
that after only 10'^ — 10^ orbits these WIMPs will undergo 
a second scattering event in the Sun. It is straightforward 
to show that this second scattering event significantly 
reduces the semi-major axis of the WIMP, removing them 
from the population of interest here (They will end up in 
the core of the Sun where they will ultimately annihilate 
with each other). 

The only way to save some of these WIMPs from this 
demise is to consider a fraction of WIMPs that have per- 
ihelions Tmin in a small range near the radius of the Sun, 
say -Rs(l ~ £) ^ ^min < Rs- Focussing on such a subpop- 
ulation of WIMPs has two advantages: (i) they traverse 
a small fraction of the mass of the Sun and therefore their 
lifetime on such grazing orbits is greatly increased, and, 
more importantly (ii) during this time, gravitational per- 
turbations due to the planets can push them into orbits 
that no longer cross the Sun. 

We first recall some concepts and notation of Hamilto- 
nian dynamics. In standard position-momenta variables 



the Hamiltonian describing the basic interaction between 
a WIMP and the Sun reads (for r = |x|), 'Hs(x,p) = 
i p2 - U{r) , where U{r) ^+GnJ p(x') d^x'/jx - x'| is 
the (spherically symmetric) Newtonian potential gener- 
ated by the mass distribution of the Sun. Note that the 
mass mx of the WIMP has been factored out of all quan- 
tities. We work with action-angle variables ("Delaunay 
variables"), traditionally denoted {L, G, H] £, g, h) [ p2[ . 
The action variables L, G, H are related to E, J and Jz, 
respectively. The angle variables (with period 2tt) corre- 
sponding to L, G, H are respectively denoted £, g, h. In 
these variables the Hamiltonian depends only on L and 
G and the general evolution equations, 
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tell us, in the case of the Hamiltonian above, that the 
action variables L, G and H are constant, while, among 
the angle variables, h is constant, but i and g evolve 
linearly in time: i = nt-\-£o, g = Co t-\-go. Here, n = 2n/P 
is the mean angular frequency of the radial motion {P ~ 
perihelion to perihelion period), and uj is the mean rate 
of advance of the perihelion. 

The crucial point to realize is the following. If we con- 
sider a WIMP orbit with a generic perihelion rmin < Rs, 
it will undergo a large perihelion precession A w ^ 27r 
per orbit, i.e. lu ^ n, because the potential U{r) within 
the Sun is modified compared to the exterior l/r po- 
tential leading to the absence of perihelion motion. In 
other words, the trajectory of the WIMP will generi- 
cally be a fast advancing rosette. This means that both 
angles £ and g are fast variables. When adding in the 
small perturbing effect due to the planets, i.e. when 
considering the total Hamiltonian, Titot = ns{L,G) + 
np{L, G, H; £, g, h; Lp, . . . ,£p, . . .) where Tip (which con- 
tains a small factor /ip = Tipianot/^s) is the planetary 
perturbation, we can work out the (first-order) secular 
effects due to the planets by averaging over the fast vari- 
ables £ and g (as well as the mean anomalies £p of the 
planets). Then the evolution equations tell us immedi- 
ately that the corresponding action variables L and G 
are secularly constant because planetary perturbations 
average out to zero (e.g. (dG/dt) = —{dn/dg)i^g = 0). 
L is essentially related to the semi-major axis a of the 
WIMP orbit, while G/L is related to the eccentricity e. 
As a result, when the rosette motion is fast, planetary 
perturbations do not induce any secular evolution in the 
semi-major axis and in the eccentricity of the WIMP or- 
bit. Such WIMPs will end up in the core of the Sun. 

A new situation arises for WIMP orbits that graze 
the Sun, because these orbits feel essentially a l/r po- 
tential due to the Sun, so that their rosette motion 
will be very slow. Consequently, the variable g will 
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be slow (compared to £), and we cannot average over 
g. We can split the total Hamiltonian in three parts 
Ti-tot — Ti-o + Ti-i + Tip where we take as unperturbed 
Hamiltonian the one corresponding to a point-like Sun 
while the perturbations are TLi = —S U (r) and Tip = 

-G^-P {j^^ - ^) ■ Here, S U{r) ^ U{r) 
GNTTT-s/r is the non-l/r part of the potential generated 
by the Sun, and Tip denotes the planetary perturbations 
involving a sum over the planets with masses Wp and 
heliocentric positions Xp. [The last term comes from 
the transformation between inertial (barycentric) coor- 
dinates and heliocentric ones.] 

Using Delaunay variables defined by the Hamiltonian 
Ho, we can derive the secular evolution of L, G, H, un- 
der the combined influence of the perturbations Tii and 
Hp. This is simply obtained by averaging the canonical 
equations over the fast angles, i.e. all the mean anoma- 
lies of the problem: i, £p. [We denote this average by 
an overbar.] By averaging the evolution equations one 
finds that (in first order) L will be secularly constant 
(i.e. a = const), while G, H, g, h slowly evolve under 
the averaged perturbed Hamiltonian 7iport(-^, G, H; g) — 
ni{L,G)+np{G,H;g;Lp). 

Because Tipert does not depend on the angle its con- 
jugate momentum H = will be secularly constant. De- 
termining the secular evolution of the remaining canon- 
ical pair {G,g), and therefore estimating the minimum 



values of the WIMP angular momenta J„ 



Gn 



for 



which planetary perturbations are strong enough to kick 
them out of the Sun is then reduced to studying the level 
curves of Hpeit{G; g). Among these level curves, there ex- 
ists a separatrix S such that initial orbits "above" S sec- 
ularly evolve into orbits with angular momenta G > Gs, 
i.e. large enough to no longer intersect the Sun. We 
find that it generically takes (for a ~ lAU) less than 
lO'^ WIMP radial periods (i.e. less than 10"^ yr) for the 
eccentricity of a WIMP above S to increase sufficiently 
to exit the Sun. Then, once G > Gs the time scale 
for the subsequent evolution of G is given by the plane- 
tary perturbations alone and is roughly l/[X]pMp(^/'^p)'^] 
longer than one orbital period, i.e. roughly 10^ yr for 
a ~ ai = 1 AU. After this time, the WIMP would, if it 
evolved only under the simplified planetary Hamiltonian 
Tip, come back again to low values of the eccentricity, 
corresponding to Sun-penetrating orbits. Under the in- 
fluence of Til , it would then again bounce back away from 
the Sun in 10^ orbits. For the scattering cross sections 
we shall be discussing below, the opacity of the small 
outer skin of the Sun is typically smaller than ^ 10^^. 
Therefore the above process could persist for hundreds 
of cycles before the WIMP gets scattered again in the 
Sun. However, it is clear that the real gravitational in- 
teraction of the WIMP with planets is much more com- 
plicated than what is described by Hp. In particular, 
the non zero eccentricities of the other planets, and the 



occurrence, once in a while, of a near collision with an in- 
ner planet will cause the elliptic elements of the WIMP 
to diffuse chaotically away from the simplified periodic 
history described above. As the very high eccentricities 
(for AU-size orbits) needed to traverse the Sun represent 
only a very small fraction of the phase space into which 
the WIMP can diffuse, on time scales of a few million 
years most of the population of WIMPs we are talking 
about will have irreversibly evolved onto trajectories that 
do not intersect the Sun for the entire age of the solar 
system. We further estimate that the long-term survival 
of such WIMPs on orbits that stay within the inner so- 
lar system is greater than 4.5 Gyr if a < ajup/2 = 2.6aE- 
[This crucial assertion is based on an approximate analyt- 
ical estimate of the lifetime of WIMPs (ll| which should 
be checked by dedicated numerical simulations.] 

We can finally estimate the density of WIMPs that 
diffuse out to solar-bound orbits and which can survive 
to the present time by simply integrating our differential 
capture rate over all trajectories J > Jmin = Gmin that 
end up out of the Sun, suitably averaged over the ini- 
tial WIMP distribution. The rate (per a = GNms/a) 
of solar capture of WIMPs that subsequently survive 
out of the Sun to stay within the inner solar system 
then depends on the A-dependent combination: gA = 
(7^ K\ , where Ja is the fraction (by mass) of ele- 
ment A in the Sun, and is the Sun-surface value of 
the capture function mentioned above. Note that the A- 
dependence is entirely contained in gA with dimensions 
[cross section] / [mass] . The total capture rate is then de- 
pendent on ^ 5^ = 5tot- 

A 

By integrating this capture rate, and making simplify- 
ing assumptions about the orbital evolution of the con- 
sidered WIMPs, we can estimate both the present space 
and velocity distributions of this new WIMP population. 
We find a local enhancement in density near the Earth, 
compared to the halo WIMP density: 
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0.21 
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where (7tot^°^ = 10^" 5tot(GeV)'^ (i.e., for masses given in 
GeV and cross sections in units of 10~^°GeV~^, with 

The a priori reasonableness of the final density en- 
hancement (5) can be seen as follows. From Eqs. (1) 
and (2) the total rate (for Jmin — 0) of WIMPs scattered 
within 1 AU is of order v^dN/da nxmsgtotv'^/vo- 
Integrating over the lifetime tg of the solar system 
and dividing by the number of galactic WIMPs present 
within lAU gives a maximum possible density enhance- 
ment (5max ~ (3/47r)(ms/a|)(i;|/'yo)is5tot ~ ^&Og[~i^°\ 
Out of this, detailed calculations of planetary perturba- 
tions of the WIMP's captured in the Sun's outer re- 
gions show lllfl that a fraction parametrically of or- 
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der (e™)"-nai/i?s)0-3(EpMp(ai/«p)') 
tracted. (Here £^=0.022 accounts for the density de- 
crease of the Sun's outer regions of interest here.) 

Perhaps a more relevant quantity, from the point of 
view of experiment, is the differential rate, dR/dQ, per 
keV per kg per day, of scattering events in a laboratory 
sample made of element A. Comparing this differential 
rate due to our new WIMP population with that due to 
the parent galactic halo population we find a ratio which 
reaches the flat maximum. 
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for energy transfers Q smaller than Qe — '^{'mx / {nT-x + 
mA))'^'rnAV^. For a target of Germanium (the present 
material of choice in several ongoing cryogenic WIMP 
detection experiments), Qe is in the keV range. If 
5tot^°^ ~ 1 this yields a 100 % increase of the differ- 
ential event rate below the value Q = Qe ^ keV which 
is typical for the energy deposit due to the new WIMP 
population, whose characteristic velocity w^; is ~ 7 times 
smaller than that of galactic halo WIMPs. 
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FIG. 1. Values of gt^^^"' as a function of WIMP nucleon 
cross section and WIMP mass for a sample range of al- 
lowed SUSY models {fj, parameter > 0) . Also shown are 
approximate experimental upper limits on WIMP nucleon 
cross section from direct detection experiments [13], assum- 
ing two different values for the local galactic halo dark matter 
density. (p — O.SGeVcm^'^ (lower) and p = O.lGeVcm^'^ (up- 
per) 

This is our central result. Its relevance depends com- 
pletely on the actual value of gtot^'^'' ■ example of 
the possibilities, in Figure 1 we display the range of gtot 



derived by sampling over the allowed Minimal Supersym- 
metric Standard Model phase space (values for SUSY pa- 
rameter /i > are chosen here), and for WIMP density 
in the range .05h~^ < px/Pciosure < 1 (where h is the 
Hubble parameter in units of lOOkms-^Mpc"^). Values 
of gtot in excess of 1 are thus possible. 

Existing detectors tend to lose sensitivity in the range 
of a few keV. However, our result provides motivation to 
consider pushing hard in this direction. The new popu- 
lation will have a strongly anisotropic velocity distribu- 
tion. Not only might this greatly help distinguish it from 
backgrounds, but a comparison of the angular anisotropy 
and annual modulation of this distribution with the cor- 
responding features of the higher energy signal from a 
halo WIMP distribution would be striking. We empha- 
size once again that if neutralinos exist in the range de- 
tectable at the next generation of detectors, this new 
WIMP population should exist at detectable levels as 
well. Finally, we note that the indirect neutrino signa- 
ture of such WIMPs that might subsequently be captured 
by the Earth, and annihilate, could be dramatic. The 
new population has a characteristic velocity that more 
closely matches the escape velocity from the Earth than 
does the background halo population. As a result, reso- 
nant capture off elements such as iron in the Earth could 
be greatly enhanced. 
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